Emergence of pygmy dipole resonances: Magic numbers and neutron skins 



Tsunenori Inakura, 1 Takashi Nakatsukasa, 1, 2 and Kazuhiro Yabana 2,1 

1 RIKEN Nishina Center, Wako, 351-0198, Japan 
2 Center for Computational Sciences, University of Tsukuba, Tsukuba 305-8571, Japan 

The pygmy dipole resonances (PDR) for even-even nuclei in 8 < Z < 40 are studied performing 
a systematic calculation of the random- phase approximation with the Skyrme functional of SkM*. 
The calculation is fully self-consistent and does not assume any symmetry in the nuclear shape of 
the ground state. In every isotopic chain, the PDR emerges by showing a peak of the El strength 
at energies less than 10 MeV. The El strength of the PDR strongly depends on the position of the 
Fermi level and shows a clear correlation with the occupation of the orbits with the orbital angular 
momenta less than 3h (I < 2). We also found a strong correlation between the isotopic dependence 
of the neutron skin thickness and the pygmy dipole strength. The fraction of the energy weighted 
strength exhausted by the PDR and the neutron skin thickness show a linear correlation with the 
universal rate of about 0.2 fm _1 . 
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Exotic nuclei show interesting new features in their re- 
sponse properties. The existence of loosely bound nucle- 
ons in the neutron skin and halo will significantly mod- 
ify the strength distributions. Among them, the low- 
energy electric dipole (El) excitation, which emerges at 
the energy well below the giant dipole resonance (GDR), 
has been most extensively studied so far. These low- 
energy peaks in the El strength distribution are often 
called the pygmy dipole resonance (PDR). The PDR 
properties have a strong impact on astrophysical phe- 
nomena. The low-energy El strengths in neutron-rich 
nuclei strongly enhance the radiative capture cross sec- 
tions for low-energy neutrons, which directly affects the 
abundance distribution in the rapid neutron capture pro- 
cess [l| . In contrast to the GDR, the PDR is sensitive to 
nuclear properties at nuclear surface and at low density. 
Thus, its property may provide us with useful constraints 
on the energy density functional, to identify the equation 
of state (EOS) of the nuclear and neutron matters. 

Another intriguing subject of interest is the relation 
between the PDR and the neutron skin thickness. The 
neutron skin thickness is known to be well correlated with 
the EOS of the neutron matter 0. Thus, if the PDR is 
correlated with the neutron skin thickness, the investiga- 
tion of the PDR may reveal the properties of the neutron 
stars, such as the proton ratio and radius. Piekarewicz 
suggested that the PDR strength is strongly correlated 
with the neutron skin thickness Q ■ This means that the 
measurement of the PDR can provide information on the 
skin thickness. In contrast, a recent correlation analysis, 
which investigates the parameter dependence of various 
properties in 208 Pb, suggests a very weak correlation be- 
tween the skin thickness and the PDR ■ These conclu- 
sions were, however, obtained from the calculations for 
specific nuclei, such as spherical Sn isotopes and 208 Pb. 
It is thus desirable to perform a systematic and fully 
self-consistent calculation for the PDR for spherical and 
deformed nuclei, to reveal the nature of the PDR and its 
relation to the neutron skin. 

The PDR were observed experimentally in neutron- 
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rich nuclei such as O isotope d, [f|, 26 Ne O 
and Sn isotopes @, but also in stable Ca 
ll and Pb isotopes [H, [3, and N = 82 isotones (TIP 
19( as well. Although the nature of the PDR, including 
its relation to the neutron skin and halo, is still elusive, 
these experiments have identified the following properties 
of the PDR: The concentrated El strengths have been 
observed around the threshold energy for the neutron 
separation. The energy-weighted strength for the PDR 
is less than 1% of the Thomas-Reiche-Kuhn (TRK) sum- 
rule value in the stable nuclei and less than 5% in the 
case of neutron-rich nuclei. 

Many microscopic studies based on the random phase 
approximation (RPA) with the Skyrme interaction [20l — 
|25| |. the Gogny interaction [2(| |2?J , and those based on 
the relativistic mean-field approach [3, [28rl3l| were car- 
ried out to investigate the PDR in the past decade. These 
microscopic calculations succeeded, to some extent, to 
reproduce dipole strength at low excitation. However, 
these microscopic calculations were applied to limited 
number of nuclei so far, and even basic questions, such 
as the condition on which the PDR appears in a nucleus, 
are not fully understood yet. 

In the present paper, we report our systematic calcu- 
lation of the El response for a wide mass region with 
8 < Z < 40 (up to A ~ 110) in the self-consistent RPA 
with the Skyrme interaction, focusing on the emergence 
and the properties of the PDR. In the present calcula- 
tion, we do not take into account the pairing correlation 
because we expect that its effect is not significant for the 
El response in light- and medium-mass nuclei [3, 25]. We 
adopt the representation of the three-dimensional (3D) 
Cartesian grids. The real-space representation has an 
advantage over other basis representations, such as the 
harmonic oscillator basis, on the treatment of the contin- 
uum, since it is well controlled by the box size and can be 
directly compared with the continuum RPA calculations. 
We adopt the 3D adaptive coordinate grids |32j within 
a sphere of a radius i?box = 15 fm. The Skyrme func- 
tional of the SkM* parameter set [33[ is used. All the 
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single-particle wave functions and potentials except for 
the Coulomb potential are assumed to vanish outside the 
sphere. The differentiation is approximated by a finite 
difference with the nine-point formula. 

The first step of the investigation is to construct the 
Hartree-Fock ground state, using the imaginary-time 
method. Then, we solve the linear-response equation for 
the El external field directly in the coordinate repre- 
sentation, avoiding explicit construction of unoccupied 
orbitals. Therefore, our calculation is free from a trun- 
cation for the unoccupied orbitals. Since the explicit 
construction of the RPA matrix in the 3D mesh rep- 
resentation is difficult, we use a new methodology, the 
finite amplitude method (FAM) [34| , in evaluation of the 
residual field Sh. The FAM allows us to construct Sh 
using the calculation of the Hamiltonian h only. The 
residual field Sh in our calculations contains all terms of 
the Skyrme interaction, i.e., the residual spin-orbit inter- 
action, the time-odd components, the residual Coulomb 
interaction, and so on. The linear-response equation is 
solved at given complex energies u) = E + using 
iterative solvers, such as the generalized conjugate resid- 
ual (GCR) method. The calculations are performed at 
energies in spacing of AE = 0.3 MeV with a fixed imagi- 
nary part 0.5 MeV, corresponding to the smearing width 
7 = 1.0 MeV. Using the obtained RPA amplitudes, we 
compute the El strength and the photoabsorption cross 
section. Details of the calculation can be found in Ref. 

The El strength of even-even nuclei are calculated up 
to zirconium isotopes from the proton to neutron drip 
lines, except for nuclei with the neutron separation en- 
ergy less than 2 MeV. This excludes the neutron drip- 
line nuclei in which extended halos develop. In addition, 
the neutron-deficient nuclei from the proton drip-line to 
N = 50 are also calculated for 40 < Z < 50. The total 
number of the calculated nuclides are 322: 40 spherical 
nuclei, 171 prolate nuclei, 56 oblate nuclei, and 55 triaxial 
nuclei. 

The calculation shows that the PDR peaks appear in 
every isotopic chain, as demonstrated in Fig. [T] for Ne 
and Ca isotopes. Increasing neutron number with a fixed 
proton number, we see that the emergence of the PDR 
suddenly takes place at N = 16 for Ne and at N = 
30 for Ca isotopes. In Ca isotopes, the PDR is rather 
distinctive and separated from the GDR peak. However, 
in the deformed neutron-rich Ne isotopes, it is not well 
separate from the low-energy tail of the GDR. Although 
it is not trivial how to define the PDR, the separate low- 
energy peaks mostly appear at energies below 10 MeV. 
Thus, in this work, the pygmy dipole strength is defined 
by the El strength at energies below 10 MeV. 

We calculate fractions of the photoabsorption cross 
section a(E) integrated up to E = 10 MeV to the in- 
tegrated total cross section. This is equal to the ratio 
TOi(PDR)/mi where mi (PDR) is the energy-weighted 
sum up to 10 MeV and mi is the energy-weighted-sum- 
rule value that is larger than the TRK value by 30 ~ 40 
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FIG. 1: (Color online) Calculated photoabsorption cross sec- 
tions in Ne and Ca isotopes. 
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FIG. 2: (Color online) Fractions of photoabsorption cross sec- 
tions mi (PDR) /mi for (a) stable and neutron-rich nuclei and 
for (b) stable and proton-rich nuclei, as functions of neutron 
number and proton number, respectively. 



%. The upper panel of Fig. [5] shows the PDR fractions 
of stable and neutron-rich nuclei as functions of neutron 
number. Here, we show isotopes with 8 < Z < 40. 
In each isotopic chain, nuclei around stable region have 
small values of the fraction less than ~1%, consistent 
with the experimental data 0, [|| [l(| • In these light nu- 
clei, a prominent PDR does not appear in stable nuclei. 
Being away from the stable region, the fractions suddenly 
increase at specific neutron numbers: N > 14, N > 28, 
and TV > 50. These "magic" numbers for the emergence 
of the PDR indicate the presence of a strong shell ef- 
fect of neutrons. The first clear indication of the kink 
structure appears at N = 14 — > 16 for O, Ne and Mg iso- 



3 



topes. For these isotopes, the strength of the PDR has a 
strong correlation with the number of neutrons occupy- 
ing the orbits with low orbital angular momenta (low-£), 
Sx/2 and d 3 / 2 (TV = 15 ~ 20). In case that the low-^ 
neutron orbits are weakly bound, they strongly expand 
in space due to the low centrifugal barrier. The kink is 
weakened by increasing the proton number and almost 
disappears for Si isotopes (Z = 14) in which TV = 16 
corresponds to the stable nucleus 30 Si. This should be 
due to the fact that these low-£ neutron orbits become 
more bound for nuclei with larger Z . The neutrons start 
filling the f?/2 orbits at TV > 20, then, the growing rate 
of the PDR strength is reduced, which is most prominent 
for Mg isotopes. 

The even more prominent kinks can be identified at 
the magic numbers of TV = 28 and TV = 50. The PDR 
fractions suddenly increase at TV — 28 — > 30 and continue 
to increase till TV — 34 where 2p shell are filled. The in- 
creasing rate of the PDR fractions depends on the proton 
number. Namely, it is the largest for small- TZ isotopes, 
such as Si and Ar isotopes, while increasing Z makes the 
rate smaller. Beyond TV = 34, the PDR fractions are 
roughly constant for 34 < TV < 50, in which the neutrons 
arc filling high-£ orbits of f$j2 and 39/2- They again show 
a sudden increase at TV = 50 — > 52, then, continue to lin- 
early grow up until 2d 5 / 2 orbits are filled at TV = 56. 
Beyond that, it is difficult to see the definite trend, since 
the ordering of the orbits and the ground-state deforma- 
tion change from nucleus to nucleus, depending on Z and 
TV. However, the careful investigation suggests that the 
occupation of s±/2 and increases the PDR fraction, 
while that of /in/2 reduces it. Thus, we may conclude 
that the spatially extended nature of the low-^ neutron 
orbits near the Fermi level plays a primary role for the 
the emergence and growth of the PDR. We have also ob- 
served that the deformation tends to increase the PDR 
strength, especially in the region TV > 56. This may be 
due to two effects; the mixture of the \ow-£ components 
in the orbits near the Fermi level and softening of the 
giant dipole resonance for the elongated direction. 

In the proton-rich side, we can see the similar behaviors 
for the isotones with 8 < TV < 50 in the lower panel of 
Fig. [2j There are kinks in the PDR fractions at Z = 14 
and 28. However, the increase as a function of the proton 
number at 28 < Z < 34 is less prominent than that as a 
function of the neutron number at 28 < TV < 34. This is 
because the Coulomb barrier prevents the low-^ proton 
orbits from spatial extension, especially in high-Z nuclei. 

Assuming the Steinwedel- Jensen model with the core 
and valence neutrons (protons), the PDR fraction is pro- 
portional to N v /N (Z v /Z) where TV„ (Z v ) is the number 
of valence neutrons (protons) [36] . This is consistent with 
our result shown in Fig. [2] for valence nucleons defined 
as those in the \ow-£ orbits. However, the model predicts 
the PDR peak energy as -Epdr/^gdr — \/N v /N, which 
disagree with our microscopic result. Figure Q] suggests 
that the PDR peak energy is rather insensitive to the 
number of nucleons. 
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FIG. 3: (Color online) Correlations between fractions of pho- 
toabsorption cross section and neutron skin thickness. The 
kink positions at TV = 28 and 50 are denoted by circles. See 
text for details. 



Next, let us examine the correlation between the PDR 
and the neutron skin thickness. The skin thickness is de- 
fined by the difference in radius between neutrons and 

protons: Ai? rms = y/ (r%) - yj (r|). We find that Ai? rms 

increases linearly with respect to the number of valence 
neutrons in the regions of TV = 14 — 20, TV = 28 — 34, 
and TV = 50 — 56, where the PDR fraction also shows 
significant linear increase as a function of the neutron 
number (Fig. [5]). Plotting the PDR fraction as a func- 
tion of Ai? rms , we then observe a linear correlation be- 
tween them. This is illustrated in Fig. |3]for isotopes with 
Z = 16 ~ 40 which show the kinks at TV = 28 and 50. 
The PDR fraction in each isotopic chain shows a linear 
correlation with the skin thickness in the regions of the 
neutron number TV = 28 — 34 and TV > 50. The positions 
of the kinks arc located at different values of Ai? rms for 
different isotopes. However, the slope is universal for all 
the isotopes; 0.18 ~ 0.20 fm -1 . It is remarkable to see 
that, even though the PDR fractions show no linear de- 
pendence on the neutron number at TV > 56 in Fig. [21 
they still keep the linear correlation with the skin thick- 
ness, until they start to decrease due to the dominant 
role of the hn/ 2 orbits at TV > 76. Despite the fact that 
the deformation and shell ordering are different and vary 
from nucleus to nucleus, the universal linear correlation 
remains valid for 50 < TV < 76. It should be noted that 
the linear correlation can be observed only for each iso- 
topic chain. Deleting the lines connecting isotopic chains 
in Fig. [3l we only see scattered points showing a weak 
correlation. 

A detailed comparison with the available experimen- 
tal data is beyond the scope of the present paper, since 
we use a smearing width of 7 = 1 MeV which prevents 
us from quantifying small low-energy El strengths which 
carry less than 1 % of the TRK sum value. So far, the 
observed PDR strength in nuclei with Z < 40 are very 
small (typically less than 1 % of the TRK value), which 
is qualitatively consistent with our calculation because 
these nuclei are in the region before the kinks (TV < 14, 
28, 50). The only exception is 26 Ne Q. The experimental 
data for 26 Ne indicates the PDR peak at E PDR ss 9 MeV 
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carrying about 5 % of the TRK value. These values cor- 
respond to our calculated values of -Epdr ~ 8 MeV with 
2.6 %. The low-energy El strength in 68 Ni is observed 
around E = 11 MeV carrying about 5 % of the TRK 
value Q. Our calculation predict a PDR peak around 
E = 10 MeV and toi(PDR)/toi indicates 4.6 % of the 
TRK value where mi (PDR) is defined by the strength 
summed up to E — 12 MeV. 

Finally, we comment on dependence of the results on 
the box size i?box of the calculation. We have selected 
several nuclei to check how the result depends on i?box- 
The discretized continuum states, which depend on the 
choice of i?box, may affect the precise peak position and 
shape of the PDR, however, the PDR strength is very lit- 
tle affected. For instance, the enlargement of the i?box to 
25 fin shifts the PDR fractions, 1.89% -> 1.94% for 26 Ne 
and 4.79% 4.65% for 60 Ca. Thus, the features of the 
PDR discussed above are valid in the present calculation. 

In summary, we carried out the systematic calcula- 
tion of the low-lying El strengths in even-even nuclei 
up to zirconium isotopes, including both spherical and 
deformed nuclei for the first time, by means of the fully 
self-consistent RPA with the Skyrme interaction utiliz- 
ing the FAM. We found the PDR in every isotopic chain. 
The PDR strength in neutron-rich nuclei shows strong 
enhancement in N > 14 for nuclei with 7 < Z < 13, in 
N > 28 for 13 < Z < 23, and in N > 50 for 23 < Z < 40. 
This suggests that the occupation of the \ow-£, such as 



s, p and d orbits, play a key role for the emergence of 
the PDR in neutron-rich nuclei. In the proton-rich side, 
the PDR strengths are hindered because of the Coulomb 
barrier to prevent those orbits from the spatial extension. 
The PDR distributions show a remarkable linear corre- 
lation with the neutron skin thickness in the regions of 
28 < N < 34 and 50 < N < 76, with a universal slope 
of 0.18 ~ 0.20 /fm. The present result is consistent with 
the linear correlation reported in Ref. [|| for neutron- 
deficient Sn isotopes in the region of 56 < N < 70, where 
the neutrons are filling the s and d orbits. Since heavy 
isotopes possess the neutron skin even in the neutron- 
deficient side, it is of great interest to perform a system- 
atic study on the PDR in heavy systems. 
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